Linear stability analysis of the Hele-Shaw cell with lifting plates 
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The first stages of finger formation in a Hele-Sfiaw cell witli lifting plates are investigated by means 
of linear stability analysis. The equation of motion for the pressure field (growth law) results to be 
that of the directional solidification problem in some unsteady state. At the beginning of lifting the 
square of the wavenumber of the dominant mode results to be proportional to the lifting rate (in 
qualitative agreement with the experimental data), to the square of the length of the cell occupied 
by the more viscous fluid, and inversely proportional to the cube of the cell gap. This dependence 
on the cell parameters is significantly different of that found in the standard cell. 



PACS number(s); 47.20.-k, 68.10.-m 
I. INTRODUCTION 

Despite of the great effort devoted lately to improve 
the understanding of Saffman- Taylor (ST) instabilities 
10,^, many related experimental facts still lack a sound 
explanation Here we are interested in the experi- 
mental observations on Hele-Shaw (HS) cells with lifting 
plates [Q. This variation with respect to the standard 
constant gap HS cell was suggested as a way to bring the 
ST problem closer to directional solidification (HJ^. In 
this experiment (see Fig. |l), instead of applying pressure 
to the less viscous fluid, the upper plate is lifted at the 
less viscous side (commonly air) at a fixed rate. It seems 
clear that the lifting of the upper plate will promote a 
pressure gradient analogous to the temperature gradient 
present in directional solidification. An interesting varia- 
tion of this experiment is the Hele-Shaw cell with a small 
gap gradient investigated by Zhao et al . 

The main results of the work reported by Ben- Jacob 
et al [Q concern the spacing of the dendrites. These 
authors observed that this spacing decreased as the lifting 
velocity was increased. They also noted a dependence on 
the initial plate spacing. On the other hand, although 
the experimental results for the HS with lifting plates 
have been analysed by means of a simplified version of 
the growth law [^,^ , a detailed analysis of this system is 
still lacking. 

In this work we present a study of growth instabilities 
in Hele-Shaw cells with lifting plates. We first derive the 
basic equation (growth law), which results to be differ- 
ent from that proposed by Ben-Jacob et al Q. In par- 
ticular our equation is not homogeneous, and thus the 
comparison with the simpler version of the directional 
solidification problem is not so evident. Then we dis- 
cuss the boundary conditions and carry out the linear 
stability analysis. Our results qualitatively explain the 
experimental data reported by Ben-Jacob et al |^ . 



II. BASIC EQUATIONS AND BOUNDARY 
CONDITIONS 

Flow in the Hele-Shaw cell is governed by the Navier- 
Stokes equation [|o|jl|] 



^ + (v.V)v = -ivP+^V2v, 
ot p p 



(1) 



where v is the speed of the fluid and P the pressure. 
p and fj, are the density and the viscosity of the fluid, 
respectively. For small Reynolds numbers and assum- 
ing that the time derivative of the fluid velocity is much 
smaller than its spatial derivatives, this equation reduces 
to 



VP = pW^v 



(2) 



Averaging over the the cell gap leads to the Poiseuille- 
Darcy equation for the mean velocity of the fluid. 
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where Vj , pj , and Pj are the velocity, viscosity and pres- 
sure field of fluid j (j=l,2); and b is the gap of the cell. 
In order to obtain an equation for the pressure field we 
need to combine Eq. (3) with mass conservation. In the 
present case the latter deserves a careful consideration. 
Let the HS cell lie in the x — y plane, the y-axis being 
the direction of motion of the fluids, and the origin of 
coordinates be at the closed (fixed) end of the cell (see 
Fig. 1^). The gap of the cell varies as 



6(y, t) = bo + y ta,n{ujt) , 



(4) 



where to is the lifting angular speed (we shall hereafter 
call a = tan^ut)). As a consequence, the mass within a 
thin column of height b changes as, Sm/5t cx 5b/dt (where 
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the density of the fluid p is assumed to be constant). 
Then, the equation which describes mass conservation 
reads 



dh 
'Ft 



-V • (6v,) 



(5) 




1. Schismatic representation of the Helle-Shaw cell 
with lifting plates, showing the parameters used in the text. 

It should be here noted that mass (and density) con- 
servation requires that, neither bubbles are formed nor 
drops of the displaced fluid are left behind during the 
lifting process. This condition, although unlikely accom- 
plished in actual experiments, it is probably unavoidable 
in analytical calculations. Eq. (5), combined with the 
Poiseuille-Darcy equation (Eq. (3)), gives the differen- 
tial equation (growth law) which governs flow in the HS 
cell with lifting plates, 
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where d expresses the derivative of a with respect to time. 
It is interesting to note that this is an inhomogeneous 
equation as opposed to the homogeneous one reported in 
Ref. [4]. The inhomogeneous term comes from the time 
dependence of the gap (Eq. (4)) and it is not present 
in the growth law for the HS with a gap gradient, as 
already found by Zhao et al Note, however, that 

the inhomogeneous term is essential in this case, as no 
fluid is injected at the open end of the cell to promote 
fluid motion (see below). We also note that including 
the spatial dependence of the gap b in Eq. (5) leads lo 
a factor 3 in the second term of Eq. (6), as opposed to 
the factor 2 reported in Ref. [4] and in agreement with 
Zhao et al 0. This change in the growth law makes 
calculations slightly more complicated. As regards the 
comparison with the directional solidification problem, 
we note that Eq. (6) is similar to that which describes 
directional solidification in some unsteady-state . 

In discussing the boundary conditions let L be the 
length of the cell and Li the length of the zone occupied 
by the more viscous fluid at t=0. On the other hand, 
and in order to carry out the linear stability analysis, we 
assume that the interface between the two fluids, instead 
of being flat, is slightly perturbed as yp — yt + Se^"^^, 
where 6 — Sqb"'* is the time-dependent amplitude of the 



perturbation (assumed to be much smaller than y^) and 
q its wavenumber. 

The first boundary condition accounts for the fact that 
the fixed end of the cell (y = 0) is closed, 



ViyiO) = 



(7) 



On the other hand, the velocities of the two fluids must 
be equal at the interface 



Vlyiyp) = V2y{Vp) 



(8) 



Finally the forces which the fluids exert on each other at 
the common interface must be equal and opposite. This 
condition can be approximated by 



Piivp) - P2{yp) 



1 



1 



RxiVjy) RziVp) 



(9) 



where a is the surface-tension (the interfacial tension as- 
sociated to the interface between the two fluids) and Rx, 
Rz , the principle radii of curvature of the interface at a 
given point. 
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(10b) 



The latter equation is related to wetting effects. Al- 
though, as already reported by other authors |]7|, |ll| , |l2| , 
we have found that these effects give a negligible contri- 
bution, we have kept this term all throughout the calcu- 
lation. 



III. LINEAR STABILITY ANALYSIS 

Once the plane interface is perturbed, the most general 
solution of Eq.(6) can be written as 



P,ix,y) = f,{y)+g,iyy''^ 



(11) 



where gjijj) is proportional to the amplitude of the per- 
turbation 6. Introducing Eq. (11) into the growth law 
(Eq. (6)) we obtain. 
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The boundary conditions for the fluid velocity at the 
closed end give 



dh 
dy 



0, 



(13a) 
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dgi 
dy 



(13b) 



Note that in order to obtain physically meaningful so- 
lutions we should also require that the perturbation be 
damped in the open end of the cell, this means that 
{dg2/dy)\y=L = 0. This assumption will be valid when- 
ever the viscosity /X2 is very small (note that, in the 
present work, wo will analyze experiments in which fluid 
2 is air). The continuity of the velocity at the interface 
between the two fluids leads to (up to first order in 6): 
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y=y, M2 V dy 



dg2 
dy 



(14a) 



• (14b) 



Finally, the continuity of forces at the interface gives 
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[9i-92{y)] 
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(15a) 
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The solution of Eq. (12a), using the boundary condi- 
tions Eqs. (13a) and (14a) can be written as. 



12/Ltjd 
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+Aj, (16) 



where Aj (j=l,2) are constants that can be determined 
from the boundary condition of Eq. (15a) (we will not 
need them here). It should be noted that, had Eq. (6) 
been homogeneous, the solution of Eq. (12a) would have 
been, fj = Aj — Bj/b{y,t)^, which, after proper use of 
the boundary conditions, gives a uniform pressure field 
and, therefore, no fluid motion, as remarked above. On 
the other hand the solution of Eq. (12b) is 



gj = j[Cjh{0+DjK^{0] 
and its derivative, 

dgj q 



dy ^ 



i = l,2, (17) 



i = l,2, (18) 



where ^ = qh{y,t)/a, and h{£,), hiO, and K^i^), K^iO 
are the flrst and the second order modifled Bessel func- 
tions. The constants Cj , Dj have to be determined from 
the boundary conditions. 



^1/2(^0)- Ai^2(Co) =0 
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K2{^i) = 0, (19c) 



(Ci - C2)/i(6) + {Di - D2)Ki{^i) = eS 



(19d) 



where ^0 = qbo/a, = q{bo + ayi)/a, = q{bo + aL)/a, 
and e is given by. 
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where the velocity of the interface between the two fluids 
is given by, 



Vi{t) = 



Of, 
12 dy 



= yi = - 



2b{yi,t) ' 



(21) 



independent of 7. Note that, due to the choice of the ori- 
gin of coordinates, this velocity is negative. Integrating 
this velocity gives an expression for the position of the 
interface (where the boundary y{t = 0) = Lj is taken): 



y^it) 



y/bl + 2Liboa - bo 



(22) 



which is the same result as if we had considered the sim- 
pler mass conservation as ^ayf + boyi — boL. 

Now, we have all the ingredients to calculate the in- 
stantaneous velocity of the interface (growth rate), which 
can be calculated from the Poiseuille-Darcy equation. 
The result is, 



viyiyp) = yp = yi + Se''^'' 



b d{fi+gie"i-) 



(23) 



12/xi dy 

Prom which we obtain the velocity of the perturbation 7, 
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where, 



^(6,Cr) 



X Ml 



-M2 



[G2(eo) + Gi(C.)] 
[G2(Co) - G2m 

[G2(a) + Gi(6)] 

[G2(a)-G2(6)] 



and, 



1,2. 



(25) 



(26) 



At the very beginning of liftting 1 << ^0 ^ ^ Cl- 

As a consequence, Gi(^i) ~ G2(Ci)- On the other hand, 
and due to the exponential behavior of the Bessel and 
Hankel functions, G2(Co) » G2(^i) » G2(^z,). Thus, 
Eq. (23b) can be approximated as 
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(27) 



The wavenumber of the dominant mode qd resuhs to 



be, 



2 _ - ^^2)v^{t) 

Qd - 
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(28) 



This result is quite similar to that obtained for the stan- 
dard Hele-Shaw cell. The only difference resides upon the 
fact that the interface velocity Vi{t) and the gap b{yi,t) 
now depend on time. As regards the cutoff wavenumber, 
we note that if the first term in the r.h.s of Eq. (27) is 
neglected (in fact it is quite small for most of the experi- 
mental configurations and conditions), the result is again 
equivalent to that of the standard cell {q^ = Sq^)- If it is 
not neglected, a minimum wavenumber, below which the 
system is stable, is also found. 
At t = 0, the dominant mode is 



,1 - fc^^.(o) 



(29) 



where Vi{0) is the velocity of the interface a.t t — 0, 



0^ 
2bn 



(30) 



Vi{0) introduces a dependence on the cell parameters of 
the wavenumber of the dominant mode, not present in 
the standard cell. In particular, qd depends on the square 
of the length occupied by the displaced (more viscous) 
fluid. It is also inversely proportional to the gap, lead- 
ing to a 6o ^ dependence of qd, as opposed to fcg ^ in the 
standard cell. These differences could be easily checked 
exp er iment ally. 



IV. NUMERICAL RESULTS AND DISCUSSION 

At the beginning of lifting (< = 0) the process is gov- 
erned by Eq. (30). In actual experiments, as the less vis- 
cous fluid is commonly air, /i2 ~ 0, and the wavelength 
of the dominant and the cutoff modes are given by ex- 
pressions identical to those for the standard cell, namely, 
Ad/^o = 7r(o-/Mi^i(0))^^^j and Ac = Xd/V^, respectively. 
In order to estimate the wavelength of the dominant 
mode we consider the experimental set up investigated 
by Zhao et al . These authors used air/glycerine (/i = 
65 mPa and (7=29.5 N/m) in a cell with a gap of 2.5 mm 
and a length of 1 m (we assume that glycerine fills the 
whole cell). Taking lu = 0.001 rad/s. Ad « 1.2 cm. This 
result decreases in an order of magnitude if, as done by 
Ben- Jacob et al Q], the lifting rate is increased in two 
orders of magnitude. We cannot precisely compare our 



results with the data obtained by the latter authors as 
they do not give important parameters such as the length 
and gap of the cell. 
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FIG. 2. Results for the velocity of the perturbation 7 
(rad/s) as a function of its wavenumber q (cm~^), for several 
times after lifting is initiated: t=Q s (continuous line), t=0.5 s 
(chain line) and t=l s (dashed line) . The results correspond 
to the cell parameters of Ref. [7], namely: air/glycerine (/i 
= 65 mPa and ct=29.5 N/m) in a cell with a gap of 2.5 mm 
and a length of 1 m (we assume that glycerine fills the whole 
cell). The lifting rate is = 0.001 rad/s. 

If instead of glycerine we consider water, at t = and 
for the same cell parameters, the wavelength of the dom- 
inant mode results to be « 0.4 cm^^, which leads to a 
wavelength of « 16 cm. This wavelength is of the order of 
the full width of standard cells and therefore one should 
not expect the formation of fingers. Further we note that 
the velocity of the perturbation results to be at least an 
order of magnitude smaller than for glycerine, reinforc- 
ing our view that the flat front would be rather stable. 
Note that this result is also obtained for the standard cell 
and that it is in agreement with the experimental obser- 
vations. The only way to increase the tendency towards 
instabilities would be to increase the velocity of the in- 
terface which in the present case can be accomplished by 
increasing L and the lifting rate and/or decreasing Bq. 

As lifting proceeds, fingers are being formed and a lin- 
ear stability analysis such as that carried out here is no 
longer strictly valid. However, some information valid to 
understand specific features of the growing process may 
be obtained from the results of a linear stability analysis, 
as already found in the case of growth in systems gov- 
erned by Poisson's equation . In Figure 2 we report our 
results for the velocity at which the perturbation propa- 
gates as a function of the perturbation wavenumber, for 
several times after lifting was initiated. In the calcula- 
tions the first term in the r.h.s. of Eq. (27) was neglected. 
The results correspond to the cell parameters given in the 
preceeding paragraph and a lifting rate of cj=0.001 rad/s. 
It is noted that the velocity of the perturbation 7 strongly 
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decreases with time. In fact the value at its maximum is 
reduced in more than a factor of 5 after lifting the cell 
for 1 second. The wavanumber at which 7 shows a max- 
imum (dominant mode) and cutoff wavenumber do also 
decrease with time (see Fig. 3). The decrease of the ve- 
locity of propagation of the perturbation suggest a lower 
tendency towards instabilities and, thus, a larger frac- 
tal dimension of the aggregates (this is compatible with 
a decreasing qd). A variation in the fractal dimension 
of the aggregates as growth (lifting) proceeds was also 
found by Roche et al Q in their simulations of the Pois- 
son equation, which they argue to be adequate for the 
HS with lifting plates (note that this is approximately 
correct if the term proportional to the first derivative of 
the pressure in Eq. (6) is either neglected or replaced 
by a constant). However they found that the simulation 
which seemed to describe more closely the experimental 
situation gave an increasing fractal dimension. 
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FIG. 3. Wavenumbers (cm~^) of the dominant (qa, con- 
tinuous line) and cutoff (qc, broken line) modes as a function 
of time (s) in the Hele-Shaw with lifting plates investigated 
in this work. The results correspond to the parameters given 
in the caption of Fig. 1. 
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Summarizing, we have presented a linear stability 
study of the HS with lifting plates. We have first de- 
rived the basic equations which result to be that of the 
directional solidification problem under some unsteady 
conditions. Despite the simplifications made to do the 
analytical work (we do not treat completely the wetting 
effects and neglect the fluid that stays attached to the 
plates), the results for the wavelength of the dominant 
mode seem to be compatible with the available experi- 
mental data. It will be worth carrying out more exper- 
iments in order to get more information on the depen- 
dence of the characteristics of the growing aggregates on 
the parameters of the cell. 
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